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Abstract 

We discuss production of two pairs of cc within a simple formalism of double-parton scattering 
(DPS). Surprisingly very large cross sections, comparable to single-parton scattering (SPS) contri- 
bution, are predicted for LHC energies. Both total inclusive cross section as a function of energy 
and differential distributions for y/s are shown. We discuss a perspective how to identify the double 
scattering contribution. 
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I. INTRODUCTION 



It is commonly believed that gluon-gluon fusion is the dominant mechanism of c or c 
production at high energies. Then in leading-order (LO) approximation the differential 
cross section for the single-parton scattering (SPS) production of heavy quark and heavy 
antiquark pair reads: 

da 1 



d yi dy 2 d 2 p t = 16^P XMX1 ^ } X29{X2 ' " } l^-^l 2 ' ( L1 ) 

where longitudinal momentum fractions can be calculated from kinematical variables of final 
quark and antiquark as: x\ = ^(exp(yi) + exp(y 2 )), x 2 = ^=(exp(— y±) + exp(— y 2 )) with 

y's being quark (antiquark) rapidities and mt being a quark (antiquark) transverse mass. 
We have limited here to gluon-gluon fusion only which is the dominant mechanism at high 
energies. The quark-antiquark annihilation plays some role only close to the kinematical 
threshold and/or large rapidities. In general, the higher-order corrections do not change 
most of observables leading to a rough renormalization of the cross section by the so-called 
K factor. 

In the present paper we wish to estimate the contribution of double-parton scatterings 
(DPS). The mechanism of double-parton scattering production of two pairs of heavy quark 
and heavy antiquark is shown in Fig. [TJ 

The double-parton scattering has been recognized and discussed already in seventies and 
eighties The activity stopped when it was realized that their contribution at those 

times available center-of-mass energies was negligible. Several estimates of the cross section 
for different processes have been presented in recent years [Tol- 18] . The theory of the double- 



parton scattering is quickly developing (see e.g. [19H26|). 

In the present analysis we wish to concentrate on the production of (cc)(cc) four-parton 
final state which has not been carefully discussed so far, but, as will be shown here, is 
particularly interesting especially in the context of experiments being carried out at LHC 
and/or high-energy atmospheric and cosmogenic neutrinos (antineutrinos). 

The double-parton scattering formalism proposed so far assumes two single-parton scat- 
terings. Then in a simple probabilistic picture the cross section for double-parton scattering 
can be written as: 

a DPS (pp -> ccccX) = -^—a SPS (pp -> ccX^ ■ a SPS (pp -> ccX 2 ). (1.2) 

This formula assumes that the two subprocesses are not correlated and do not interfere. 
At low energies one has to include parton momentum conservation i.e. extra limitations: 
X\ +£3 < 1 and x 2 + x 4 < 1, where X\ and X3 are longitudinal momentum fractions of gluons 
emitted from one proton and x 2 and X4 their counterpairs for gluons emitted from the second 
proton. The "second" emission must take into account that some momentum was used up 
in the first parton collision. This effect is important at large quark or antiquark rapidities. 
Experimental data [27] provide an estimate of a e ff in the denominator of formula (jl.2p . In 
our analysis we take a rather conservative value a e ff = 15 mb. 

The simple formula ( II. 2p can be generalized to include differential distributions. Again 
in leading-order approximation differential distribution can be written as 

da 1 da da 

(1.3) 



dy 1 dy 2 d 2 p lt dy 3 dy 4 d 2 p2t 2a eff dy 1 dy 2 d 2 p lt dy 3 dy 4 d 2 p 2t 
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FIG. 1: Mechanism of cccc production via double-parton scattering. 



which by construction reproduces formula for integrated cross section (11. 2ft . This cross sec- 
tion is formally differential in 8 dimensions but can be easily reduced to 7 dimensions noting 
that physics of unpolarized scattering cannot depend on azimuthal angle of the pair or on 
azimuthal angle of one of the produced c (c) quark (antiquark). The differential distributions 
for each single scattering step can be written in terms of collinear gluon distributions with 
longitudinal momentum fractions x±, x 2 , x 3 and x 4 expressed in terms of rapidities yi, y 2 , y 3 , 
?/4 and transverse momenta of quark (or antiquark) for each step (in the LO approximation 
identical for quark and antiquark). 

A more general formula for the cross section can be written formally in terms of double- 
parton distributions, e.g 
at high energies: 



F gg , F qq , etc. In the case of heavy quark (antiquark) production 



da 



DPS 



1 



2(T eff 

dcr gg ^cc.^S^\i *^*x? f^x^^^QQ ^cc(^2; *^2' A^2/ dx \ dx% doc ^ dx 2 • 



;i-4) 



It is physically motivated to write the double-parton distributions rather in the impact 
parameter space F gg (xi,X2,b) = g{x\)g{x2)F{b) 1 where g are usual conventional parton 
distributions and F(b) is an overlap of the matter distribution in the transverse plane where 
b is a distance between both gluons in the transverse plane 28) . The effective cross section 



in ( II. 2p is then l/cr e // = J d 2 bF 2 (b) and in this approximation is energy independent. 

The double-parton distributions in Eq. (ll.4|) are generally unknown. Usually one assumes 
a factorized form and expresses them via standard distributions for SPS. Even if factorization 
is valid at some scale, QCD evolution may lead to a factorization bre aking. Evolution is 

for 
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i.e. 



known only in the case when the scale of both scatterings is the same [19|, l 
heavy object, like double gauge boson production. For double cc production this is not the 
case and was not discussed so far in the literature. In the present preliminary study we shall 
therefore apply the commonly used in the literature factorized model. A refinement will be 
done elsewhere. In explicit calculations presented below we use leading order collinear gluon 
distributions (GRV94 CTEQ6 0, GJR08 H3, MSTW08 H). 
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II. RESULTS 



In Fig. [2] we compare cross sections for the single and double-parton scattering as a 
function of proton-proton center-of-mass energy. At low energies the conventional single- 
parton scattering dominates. For reference we show the proton-proton total cross section as a 
function of energy as parametrizes in Ref . 33[ • At low energy the cc or cccc cross sections are 
much smaller than the total cross section. At higher energies the contributions dangerously 
approach the expected total cross section 1 . This shows that inclusion of unitarity effect 
and/or saturation of parton distributions may be necessary. The effect of saturation in cc 
production has been included e.g. in Ref. 34| but not checked versus experimental data. 
Presence of double-parton scattering changes the situation. The double-parton scattering 
is therefore potentially very important ingredient in the context of high energy neutrino 



production in the atmosphere [34H36I] or of cosmogenic origin [37J. We leave this rather 
difficult issue for future studies where the LHC charm data must be included. At LHC 
energies the cross section for both terms become comparable 2 . This is a completely new 
situation when the double-parton scattering gives a huge contribution to inclusive charm 
production. 




FIG. 2: Total LO cross section for single-parton and double-parton scattering as a function of 
center-of-mass energy (left panel) and uncertainties due to the choice of (factorization, renormal- 
ization) scales (right panel). We show in addition a parametrization of the total cross section in 
the left panel. Cross section for DPS should be multiplied in addition by a factor 2 in the case 
when all c (c) are counted. 

In Figs. El HI we present single c (c) distributions. Within approximations made in this 
paper the distributions are identical in shape to single-parton scattering distributions. This 
means that double-scattering contribution produces naturally an extra center-of-mass en- 
ergy dependent K factor to be contrasted with approximately energy-independent i^-factor 
due to next-to-leading order corrections. One can see a strong dependence on the factor- 
ization and renormalization scales which produces almost order-of-magnitude uncertainties 



1 New experiments at LHC will provide new input for parametrizations of the total cross section. 

2 If inclusive cross section for cor c was shown the cross section should be multiplied by a factor of two - 
two c or two c in each event. 
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and precludes a more precise estimation. A better estimate could be done when LHC charm 
data are published and the theoretical distributions are somewhat adjusted to experimental 
data. 
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FIG. 3: Distribution in rapidity (left panel) and transverse momentum (right panel) of c or c 
quarks at y/s = 7 TeV. Cross section for DPS should be multiplied in addition by a factor 2 in the 
case when all c (c) are counted. 




FIG. 4: Uncertainties related to renormalization and factorization scales choice for distributions 
in rapidity (left panel) and transverse momentum (right panel) of c or c quarks at y/s = 7 TeV. 
Cross section for DPS should be multiplied in addition by a factor 2 in the case when all c (c) are 
counted. 

So far we have discussed only single particle spectra of c or c (rapidity, transverse mo- 
mentum distributions) which due to scale dependence do not provide a clear test of the 
existence of double-parton scattering contributions. A more stringent test could be per- 
formed by studying correlation observables. In particular, correlations between c and c are 
very interesting even without double-parton scattering terms [38[ . In Fig. we show dis- 
tribution in the difference of c and c rapidities Udiff = Vc — Vc (left panel) as well as in 
the cc invariant mass M cg (right panel). We show both terms: when cc are emitted in the 
same parton scattering [c\Ci or C3C4) and when they are emitted from different parton scat- 
terings (C1C4 or C2C3). In the latter case we observe a long tail for large rapidity difference 
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as well as at large invariant masses of cc. Such distributions cannot be directly measured 
for cc but could be measured for mesons (rapidity difference up to 5 for the main ATLAS 
or CMS detector) or electron-positron or The ALICE forward muon spectrometer 

covers the pseudorapidity interval 2.5 < rj < 4 which when combined with the cen- 



tral detector means pseudorapidities differences up to 5. This is expected to be a region 
of phase space where double-parton scattering contribution would most probably dominate 
over single-parton scattering contribution. This will be a topic of a forthcoming analysis. 
Next-to-leading order corrections are not expected to give major contribution at large pseu- 
dorapidity differences or/and large invariant masses of fi + n~ but this must be verified in 
the future. The CMS detector is devoted especially to measurements of muons. The lower 
transverse momentum threshold is however rather high, the smallest being about 1.5 GeV 
at T) — ± (2 - 2.4) which may be interesting for double-parton scattering searches. This 
requires special dedicated Monte Carlo studies. 



: DPS pp->ccccX 

FLO: g g c c 
: CTEQ6 LO 



s = 7 TeV 



c,c 4 or c 2 c 3 
c,c, or c,c. 




> 

CD 

i5 



"D 



-i 1 1 1 1 1 r— | -r— i 

DPS pp->ccccX 

' LO: g g -> c c 



s = 7TeV : 



c,c 4 or c 2 c 3 n 



c,c 2 or c 3 c 4 - 



'diff 




FIG. 5: Distribution in rapidity difference (left panel) and in invariant mass of the cc pair (right 
panel) at y/s = 7 TeV. 



Finally in Fig. [6] we present distribution in the transverse momentum of the cc pair |p± c g|, 
where p±cc — P±c + Rls which is a Dirac delta function in the leading-order approximation. 
In contrast, double-parton scattering mechanism provide a broad distribution extending 
to large transverse momenta. Next-to-leading order corrections obviously destroy the (5-like 
leading-order correlation. We believe that similar distributions for DD or/and e + e~ or fi + fi~ 
pairs would be a useful observables to identify the DPS contributions but this requires real 
Monte Carlo simultions including actual limitations of experimental apparatus. Correlations 
between outgoing nonphotonic electrons has been studied at much lower RHIC energy in 
Ref. (!}. 

Production of two cc pairs in the leading order approximation is only a first step in trying 
to identify DPS contribution. In the next step we plan next-to-leading order calculation of 
the same process. Inclusion of hadronization and/or semileptonic decays would be very 
useful in planning experimental searches. 
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FIG. 6: Distribution in transverse momentum of cc pairs from the same parton scattering and 
from different parton scatterings at y/s = 7 TeV. 
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